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Abstract 

Many studies have been carried out in order to increase the search efficiency of con- 
straint satisfaction problems; among them, some make use of structural properties of the 
constraint network; others take into account semantic properties of the constraints, gener- 
ally assuming that all the constraints possess the given property. In this paper, we propose 
a new decomposition method benefiting from both semantic properties of functional con- 
straints (not bijective constraints) and structural properties of the network; furthermore, 
not all the constraints need to be functional. We show that under some conditions, the 
existence of solutions can be guaranteed. We first characterize a particular subset of the 
variables, which we name a root set. We then introduce pivot consistency, a new local 
consistency which is a weak form of path consistency and can be achieved in 0(n 2 d 2 ) com- 
plexity (instead of 0(n 3 d 3 ) for path consistency), and we present associated properties; in 
particular, we show that any consistent instantiation of the root set can be linearly ex- 
tended to a solution, which leads to the presentation of the aforementioned new method 
for solving by decomposing functional CSPs. 

1. Introduction 

Improving search efficiency is one of the main goals of artificial intelligence. Efficiency is 
particularly important for researchers interested in constraint satisfaction problems (CSPs), 
also called constraint networks. In order to achieve this aim, various approaches have 
been tried, such as filtering techniques (mainly arc and path consistencies) (Montanari, 
1974; Mackworth, 1977; Mohr & Henderson, 1986; Bessiere, 1994) or improvements to the 
backtrack process (Haralick & Elliot, 1980; Dechter & Pearl, 1988; Ginsberg, 1993; Prosser, 
1993). Other work concerned the characterization of classes of polynomial problems, based 
on the size of their domains (Dechter, 1992) or on the structure of the constraint network 
(Freuder, 1978), leading to the presentation of decomposition methods such as the cycle- 
cutset (Dechter, 1990) or the tree clustering (Dechter & Pearl, 1989) methods. A more 
recent approach consists of taking into account semantic properties of the constraints (as 
opposed to structural or topological properties of the network) to achieve arc consistency 
efficiently for specific classes of constraints (Van Hentenryck, Deville, & Teng, 1992; Mohr 
&: Masini, 1988), or to characterize some tractable classes of problems (van Beek, 1992; van 
Beek & Dechter, 1994; Kirousis, 1993; Cooper, Cohen, & Jeavons, 1994). We present these 
results further in this paper. 

Some frequently encountered constraints are functional constraints, for instance in pep- 
tide synthesis (Janssen, Jegou, Nouguier, Vilarem, &: Castro, 1990) or in Constraint Logic 
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Programming (Van Hentenryck et al., 1992). Many constraint-based systems and languages 
use functional constraints, including THINGLAB (Borning, 1981), GARNET (Myers, Giuse, Sz 
Vander Zanden, 1992) and kaleidoscope (Freeman-Benson & Borning, 1992). 

In this paper, we study functional constraints: the relations they represent are partial 
functions (notice those are not the same as bijective constraints: see sections 2.1 and 2.2). 
More precisely, we study the properties given to them by a new local consistency, pivot 
consistency. We notably show that, under some conditions, this local consistency guar- 
antees the existence of solutions and makes them easier to find: first, we characterize a 
subset of the variables, called a root set and denoted 7Z; we then show that if the network 
is pivot consistent, then any consistent instantiation of 7Z can be linearly extended to a so- 
lution, provided the instantiation order is 1Z- compatible (i.e., it possesses some topological 
properties associated with 7Z and the functional constraints). 

We then introduce a new method for solving any functional CSP by decomposing it 
into two subproblems: first, finding a consistent instantiation of the root set, and second, 
extending this partial instantiation to a solution (in a backtrack- free manner). 

The main difficulty is clearly to find a consistent instantiation of the root set, which 
remains exponential in its size. The method we present here is therefore all the more efficient 
since the root set is of small size. 

Another aspect we wish to point out is that, unlike most of the work dealing with 
constraints which possess some specific properties, this method not only applies when all 
the constraints possess the given property (or properties), but also when only some of them 
do. 

This paper is organized as follows: section 2 first introduces the basic definitions and 
notations regarding CSPs which will be used in the following. We then define functional 
constraints and discuss previous work and its relations with this paper. The last part of the 
section presents notions associated with functional constraints, which are mainly found in 
graph theory. We then introduce pivot consistency and an 0(n 2 d 2 ) algorithm allowing to 
achieve the latter in section 3, before presenting in section 4 the properties supplied to a CSP 
by this consistency. We finally propose a method based on those properties for solving CSPs 
composed of functional constraints (with or without additional non- functional constraints). 

2. Preliminaries 

Definition 2.1 A binary CSP V is defined by (X, D,C, R), where 

• X = {X\, X2, . . . , Xi, . . . , x n } is the set of its n variables. 

• D = {di, . . . , Dj, . . . , D n } is the set of its n domains, where each Di = {aj, 6j, . . .} is 
the set of the possible values for variable Xi; the size of the largest domain of D is 
called d. 

• C is the set of its e constraints {ci, . . . , c&, . . . , c e }; for the sake of simplicity, a 
constraint c\ = {Xj,Xj} involving variables Xj and Xj is denoted c\j. 

• R is the set of the e relations associated with the constraints, where the relation Rij is 
the subset of the Cartesian product Di x Dj specifying the pairs of mutually compatible 
values. 
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A graph Q = (X, C) called the constraint graph can be associated with V, where the vertices 
and the edges (or arcs) respectively represent the variables and the constraints. Another 
graph, representing the relations of the CSP, can also be defined: the consistency graph. 
Figure 1 represents the constraint and consistency graphs of the following example: 

A travel agency offers five destinations to its customers: Paris, London, Wash- 
ington, New- York and Madrid, and employs three guides to accompany them: 
Alice, Bob and Chris. These guides provide the information below: 

• Alice wishes to go to Paris or New-York, but she only speaks French. 

• Dob speaks English and French, and does not want to go to Madrid (but 
accepts any other city). 

• Chris only speaks Spanish, and refuses to go to any city but New- York. 

The manager would like to find all the possibilities left to him, i.e., which guide 
he can send to which city, with which currency, knowing that each guide must 
of course speak the language of the country he (or she) visits, and the currency 
must be that of the country. 

This problem can be encoded by the CSP V ex , composed of five variables and five constraints: 

• The set of the five variables is 

X = {GUIDES, CITIES, COUNTRIES, CURRENCIES, LANGUAGES} 

• Their five domains are: 

Dguides= {Alice (A), Bob (B), Chris (C)} 

d cities = {Paris, London, Washington, New-York, Madrid] 

d countries = {France, GB, USA, Spain) 

d currencies = {FrF, £, $, Pes} 

T>LANGUAGES = {French, English, Spanish} 

• The set of the five constraints is: 

C = {C GUIDES-CITIES, C CITIES-COUNTRIES, C GUIDES-LANGUAGES, 
CCO UNTRIES- C URRENCIES , C CO UNTRIES-LA NGUA GEs} 

where these constraints respectively represent the cities the guides wish (or accept) 
to visit, the countries the cities belong to, the languages spoken by the guides, the 
currencies used in the countries, and the official languages of the countries. Notice 
there is no explicit constraint between the currencies and the guides: this constraint 
will be induced by the cities (and thus the countries) visited by the guides on the one 
hand, and by the currencies of the countries on the other hand. 

• The five associated relations are: 

R GUIDES-CITIES = {(Alice, Paris), (Alice, New- York), (Bob, Paris), (Bob, London), 
(Bob, Washington), (Bob, New- York), (Chris, New-York)} 

Rguides-languages = {(Alice,French), (Bob,French), (Bob, English), 
( Chris, Spanish )} 
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Figure 1: Constraint graph and Consistency graph of V ( 
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Rcities-COUNTRIES = {(Paris, France), (London, GB), (Washington,USA), 
(New- York, USA ), (Madrid, Spain )} 

Rcountries-currencies = { (France, FrF), (GB,£), (USA,$), (Spain,Pes)} 

R-co UN tries-lang UA ges = {(France, French), (GB, English), (USA, English), 
( Spain, Spanish)} 

If there is no constraint between variables Xi and Xj (Cij C), all pairs of values are 
therefore allowed: Rij = x Dj. Rij is then called a universal relation. 

The notion of support is related to the pairs that belong or not to the relations: we say 
that a value a& 6 is a support of (or supports) a value aj 6 Dj for the constraint Cj& if 
and only if (ai,at) 6 Rik- Also, the value 6 is a support of the pair (a,i,aj) 6 if 
and only if both supports ai for the constraint C\k and aj for the constraint Cjk- 

Y stands for any subset of X, and Y^ represents the set of the first k variables of X, 
w.r.t. an order explained in the context. The default order will be Xi,X2, ■ ■ ■ ,x n _i,x„. 
Also, the order for V ex will be GUIDES, CITIES, COUNTRIES, CURRENCIES, LANGUAGES. 

Given Y t = {x i, . . . , x t } a subset of X, the assignment to each variable Xj of Y t 
of a value of Dj is called a (partial) instantiation of Yt. Such an instantiation It = 

. . . , aj, . . . ,a,j, . . . ,at) is consistent if and only if every constraint included in Y t is satis- 
fied; formally, Mi,j < t s.t. c\j 6 C, (a,i,aj) 6 R^j. For instance, neither (Alice, London) nor 
(Bob, Paris, Spain) are consistent, whereas (Alice, Paris) and (Chris, New-York, Peseta, 
Spanish) both are. A solution is therefore defined as a consistent instantiation of X. A 
solution to X is thus (Alice, Paris, France, FrF, French). 
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Figure 2: Two examples of functional constraints 



2.1 Functional Constraints: Definition 

We first define what we call functional constraints (Figure 2 presents two examples of such 
constraints), before discussing related work. We then recall some graph theory notions 
before characterizing a property necessary for the method we present in this paper. 

Definition 2.2 Given two variables Xj and x^, we denote Xj — > iff for all ai in Di, 
there is at most one a& in such that (aj,a&) 6 R^. We then denote a& = fikip-i) (or 
au = fi^ki^i)) if this value exists, s otherwise. 

A constraint Cik is functional iff Xj — > x& or x^ — > Xj. In the following, we will call Xj 
the origin and x& the target of the functional constraint Xj — > x&. 

From this definition, we can deduce a partition of the constraint set of the CSP: on the 
one hand C/, the set of the functional constraints, on the other hand C„, the others (the 
non functional constraints). A CSP is said to be functional if it contains some functional 
constraints (Cf ^ 0), and strictly functional if it only contains such constraints: C = Cf 



In the field of constraint satisfaction problems, one of the most important objectives is 
to reduce computation time. Many studies have been carried out to achieve this (see for 
instance the references in the introduction of this paper). They permit to reduce, sometimes 
drastically, the time required to solve CSPs. However, an obstacle limits their efficiency: 
they are general methods, and, as such, they do not take into account the specificities of 
the problems to be solved, notably w.r.t. semantics of the constraints. This deficiency was 
overcome by various (and, for most of them, recent) results. They exhibit some classes of 
constraints which possess particular properties. These results fall into two classes: on the 
one hand, those proposing algorithms specifically adapted to some classes of constraints, 
more efficient than general algorithms, and on the other hand, those characterizing classes 
of "tractable" problems. In this section as in the rest of this paper, we will suppose that, 
unless otherwise stated, the networks are binary. 



(C o = 0)- 



2.2 Related Work 
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The first class contains those techniques which make use of the semantics of the con- 
straints to propose specific processings, namely arc consistency filterings. 

• Mohr and Masini (1988) propose to modify AC-4 to improve its efficiency on equality, 
inequality and disequality constraints. They reduce time complexity of AC-4 from 
0{ed 2 ) in general case, to O(ed) for those specific constraints. 

• AC-5, a new arc-consistency algorithm, is presented by Van Hentenryck et al. (1992). 
Actually, AC-5 is not an algorithm: it should rather be considered as a model with 
several possible specializations, according to the classes of constraints to process {e.g., 
functional 1 , anti-functional and monotonic constraints — notice that these constraints 
respectively correspond to Mohr and Masini's equality, disequality and inequality con- 
straints). AC-5 achieves arc-consistency in O(ed) for those constraints, which is con- 
sistent with the former result. 2 

The second class contains those techniques which identify constraints for which some 
local consistency is sufficient to guarantee global consistency, or propose polynomial algo- 
rithms to solve some classes of problems. 

• van Beek (1992) shows that path consistency induces global consistency if the con- 
straints are row convex. 3 

• van Beek and Dechter (1994) show that if the constraints are m-tight 4 , then strong 
(m + 2)-consistency induces global consistency. 5 

• A tractable class of binary constraints is identified by Cooper et al. (1994), 0/1/all 
constraints. If for any pair of variables Xj, Xj, every value of r>j is supported by either 
0, 1 or all values of Dj, then there exists an algorithm, called ZOA (Zero/One/All), 
which either finds a solution, or informs that no solution exists in polynomial time 
(namely 0(ed(n + d))). It is also shown that any other class of constraints closed 
under permutation and restriction of the domains is NP-complete. 

• Kirousis (1993) presents a polynomial algorithm to solve CSPs with implicational con- 
straints (which, for binary constraints, are the same as 0/1/all constraints 6 ). 

1. Which should in fact be called bijective. 

2. It is also shown that, in the field of constraint logic programming over finite domains, for a restricted class 
of basic constraints, node- and arc-consistency algorithms provide a decision procedure whose complexity 
is the same as the complexity of AC-5, namely O(ed). 

3. A binary constraint Cij is row convex iff for any at € Di, all its supports are consecutive in Dj, provided 
there exists an (implicit) ordering of the domains. 

4. A binary constraint Cij is m-tight iff every value at € D{ is supported by either at most m values in 
Dj, or all values of Dj, and, conversely, every value a,- € Dj either has at most m supports in Di or is 
supported by all values of Di. 

5. They furthermore show how this result applies to non-binary networks: for any CSP whose constraints 
have arity r or less and are at most m-tight (an r-ary constraint is m-tight iff all its binary projections 
are m-tight), if the CSP is strongly ((m + l)(r — 1) + 1) -consistent, then it is globally consistent. 

6. For the sake of simplicity we will not present non-binary implicational constraints (for which an 0(d 2 r 2 ne) 
algorithm is presented), nor a parallel algorithm, also studied in his paper. 
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A first comment about those results: a class of binary constraints is a subclass of all those 
presented above, namely bijective constraints: they are 1-tight, 0/l(/all), and row convex. 
Bijective constraints are known to be tractable. Notice that constraints called functional by 
van Beek (1992) and by Van Hentenryck et al. (1992) are actually bijective constraints 7 . On 
the other hand, a functional constraint Xj — > Xj is neither necessarily 1-tight, nor 0/1/all 
(see for example constraint w = v 2 in Figure 2), but Rij remains row convex (even though 
Rji may not be). But achieving path consistency in a network whose constraints are all 
functional may create non functional constraints (or transform some functional constraint 
into non functional ones), van Beek's result does not apply anymore (since there are non 
row convex constraints). Functional constraints, unlike bijective constraints, are intractable 
in general case. 

The second comment is more important: all of those results assume that all of the 
constraints of the network belong to a given class (e.g., m-tight, 0/1/all, row convex. . . ), 
and only apply in this case. The difference, and, in our opinion, an interesting improvement, 
is that the results we present in this paper also apply when only some of the constraints 
are functional. 

In a previous paper (David, 1993), it is shown that in a path consistent CSP, any 
consistent instantiation can be extended to a solution if, for every non-instantiated variable, 
there exists a sequence of functional or bijective constraints from an (any) instantiated 
variable to this non-instantiated variable. A subset of the variables is thus defined, such 
that any of its consistent instantiation can be linearly extended to a solution. Actually, path 
consistency is not necessary: a weaker local consistency, pivot consistency, is sufficient, as 
we now propose to show. 

2.3 Functional Constraints and Directed Graphs 

If a CSP V is functional, its constraint graph can be divided into two subsets. Cf contains the 
(directed) arcs representing the functional constraints, where an arc is directed from Xi to 
Xk iff Xj — > Xfc, and C contains the edges (i.e., undirected), for non functional constraints. 
We denote Q = (X, Cf U C ) this graph, and Qf = (X, Cf) its directed subgraph. 

Before going further, we need to recall some graph theory notions: the descendant of a 
vertex (variable), and the root of a directed graph: 

Definition 2.3 (Berge, 1970) A vertex x^ is a descendant of Xj in a directed graph 
(X, Cf) iff there exists in Cf a path (Xi,Xj 1 , . . . ,Xj , . . . ,Xj p ,x^) such that X{ — > Xj 1 , . . . , 

X U x k+v ■ ■ > X 3 P -> x k, Vg € 1-P - 1. 

A vertex x r is a root of a directed graph (X,Cf) iff any other vertex is a descendant 
from x r . 

By extension, we define a root set of a directed graph. It is easy to find out that a root is a 
root set which contains a single vertex. 

Definition 2.4 A subset TZ of X is a root set of a directed graph (X,Cf) iff any vertex of 
X — TZ is a descendant from an element of 1Z. 1Z is called minimal iff there does not exist 
a root set TZ' such that 1Z' C 1Z. 1Z is called minimum iff there does not exist a root set TV 
such that \TZ'\ < \TZ\. 

7. A constraint Cy is bijective iff it is functional from Xi to Xj (xi —¥ Xj) and from Xj to Xi (xj —¥ Xj). 
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Q = {X,C S UC ) Q f = {X,C f ) 

Figure 3: Subgraph and root set of V ex 

GUIDES, CITIES, COUNTRIES, CURRENCIES, LANGUAGES 
or CITIES, GUIDES, COUNTRIES, LANGUAGES, CURRENCIES 
are both 7£-compatible. 

GUIDES, CITIES, LANGUAGES, COUNTRIES, CURRENCIES 
is prefixed by 1Z, but it is not 7£-compatible: LANGUAGES 
has no ancestor Xi in this ordering such that Xi — > LAN- 
GUAGES. 

GUIDES, LANGUAGES, COUNTRIES, CURRENCIES, CITIES 
is not 7^-compatible: it is not prefixed by 1Z. 

Figure 4: 7^-compatible and non 7^-compatible orderings of V ex 



CURRENCIES LANGUAGES 



COUNTRIE. 




Figure 3 illustrates the notions of directed subgraph and of root set (here, a minimum root 
set) applied to our example V ex . Notice that the arc is only directed from COUNTRIES to 
CURRENCIES, and not from CURRENCIES to COUNTRIES though we have both COUNTRIES 

CURRENCIES and CURRENCIES COUNTRIES. This is only for the sake of simplicity: 
this information will not be useful in the following. 

With a root set 7Z, we associate 7Z- compatibility, a property that can be verified by an 
ordering of X: 

Definition 2.5 An ordering of X X\,Xi, . . . ,x n is called 7^-compatible iff: 

1. V« < \K\, Xi 6 1Z, and 

2. \/k > \1Z\, 3j < k such that Xj — > x^. 

In other words, the first variables in this ordering are those which belong to 1Z (we then say 
that this ordering is prefixed by 7Z), and any other variable x^ has at least one ancestor Xj 
in the ordering (xj is ordered before x^) such that Xj — > x^ (see Figure 4). 

3. Pivot Consistency 

First of all, why did we call it pivot consistency? Pivot consistency depends on the assign- 
ment ordering 8 . On every step of the assignment, two elements have to be distinguished: 
firstly, the set of the formerly instantiated variables (the instantiated set), and secondly the 

8. Similarly to directional path consistency and adaptative consistency (Dechter & Pearl, 1988). 
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next variable to be instantiated. One can consider pivot consistency as a property between 
each instantiated set and the next variable to be instantiated. For each of these sets, a 
constraint has a particular role to play: the consistency checks turn on it. This constraint 
is thus called the pivot of the set. 

3.1 Presentation 

Further in this section, we will introduce the definition of pivot consistency in 3 steps: 
consistency between 3 variables, then pivot of a subset of X, and finally pivot consistency 
of the CSP. 

Definition 3.1 Let Xi,Xj,x k G X such that C'i k ,Cj k £ C. Ci k and Cj k are x k -compatible 
iff any tuple of Rij has at least one support in D k for Ci k and Cj k . Formally, V(aj,aj) € Rij, 
3a k G D k s.t. (ai,a k ) G Ri k and (aj,a k ) G Rj k 

Four comments arise from this definition: 

1. The existence of a constraint between Xj and Xj is not compulsory: the relation R^ 
may be universal. 

2. Ci k and Cj k are not necessarily different. If not, c^ k must be x k - compatible with itself; 
in that case we assimilate the pair (ai,ai) of the relation Ra to the value ai of the 
domain D{. 

3. This definition may be seen as a local version of strong 3-consistency (2 and 3 consis- 
tencies): 

• 2-consistency: it is due to the remark above: any value in Dj must have a support 
in D k . 

• 3-consistency: any consistent instantiation of {Xi,Xj} may be extended to a 
third variable, here x k . 

4. Knowing that path consistency is equivalent to 3-consistency, we can deduce that path 
consistency can be rewritten in terms of x^-compatibility: 

A CSP is path consistent iff for all Xi,Xj,x k pairwise distinct, {Xj, x k } and {xj, x k } 
are both x k - compatible. 

We already said that some constraints called the pivots have a particular part to play: we 
now define them: 

Definition 3.2 Let Y C X and c ik G C s.t. Xi — > x k , x k 6 X — Y and x.- b 6 7. Xi — > x k 
is a pivot of Y iffVxj G Y s.t. Cj k 6 C, c\ k and Cj k are x k - compatible. 

In other words, given any proper subset Y of X , a functional constraint Xi — > x k "coming 
out" from Y (Xj 6 Y and x k G X — Y) is a pivot of Y if and only if for any consistent 
instantiation (a,i,a,j) of Xj and of any other variable Xj in Y, there exists at least one 
value a k in D k such that (ai,aj) can be extended to a consistent instantiation (ai,aj,a k ) of 
{Xi, Xj, x k }. As for functional constraints, we call Xj the origin and x k the target of the 
pivot Xi — > x k . 
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n = {x 1 ,x 2 } 

The pivot of {xi, x 2 } is x 2 -> x 3 

This CSP is pivot consistent: the 3 pairs of R\ 2 , that 
is (01,02), (61,02) and (61,02) all have a support in 
D3 (respectively 03, 03 and 63). 

On the other hand, it is neither path consistent, since 
(b 2 , 63) has no support in D\, nor arc consistent, since 
62 has no support in D\. 

Figure 5: A pivot consistent CSP 



TZ = {x 1 ,x 2 } 

The pivot candidate of {x\,x 2 } is x 2 — > x$. 

This CSP is not pivot consistent: (ai, b 2 ) has no sup- 
port in D3. 



Figure 6: A non pivot consistent CSP 

Now that the pivot of a subset Y of X has been defined, we can introduce the notion of 
pivot consistency of a CSP: 

Definition 3.3 Given a CSP V = (X, D, C,R), a root set 7Z of (X, C) and an 1Z- compatible 
assignment ordering X\, x 2, . . . , x n , we say V is pivot consistent w.r.t. this ordering iff 
Mk > r = \7Z\, 3h < k s.t. x^ — >■ x^ is a pivot of Y^-i = X\, X2, . . . , ^fe-i- 

Informally, every variable which does not belong to 1Z is the target of a pivot whose origin 
is before it in the assignment ordering. Figures 5 and 6 show a pivot consistent CSP and a 
non pivot consistent CSP. 

Pivot consistency of a CSP relies on the existence of a set of functional constraints 
possessing particular properties: the pivots (Definition 3.3 above). Thus, a minimum set of 
functional constraints can be characterized, so that the CSP is pivot consistent. This is the 
purpose of the three conditions stated in Definition 3.4: 

• condition 1 ensures there exists a pivot for each Yfc_i, k > r: the network is therefore 
pivot consistent; 

• conditions 2 and 3 leave unnecessary constraints out: only one pivot needs to "target" 
each Xk G X — TZ, and no pivot is required inside 1Z. 
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Definition 3.4 Given a CSP V = (X,D,C = Cj U C ,R), a root set TZ, an TZ- compatible 
assignment ordering X\, X2, ■ ■ ■ , x n , and a set of functional constraints P C Cf, if the 
following three conditions are satisfied 

1. Vx k E X — TZ, 3x h -» x k E P s.t. h < k and x h — > x k is pivot of Y k _ x 
(any variable of X — TZ is the target of a pivot) 

2. Vx fc E X - TZ, {x h ->x t eP and Xj -t x k E P} =>■ h = j 
(any variable of X — TZ is the target of at most one pivot) 

3. MXj E TZ, fix, -> Xj E P 

(no variable of TZ is the target of a pivot) 

then P is called a pivot set of the CSP, and V is pivot consistent. 

We have introduced a new local consistency, pivot consistency. However, as for any local 
consistency, a CSP generally does not satisfy it. The problem has then to be filtered in order 
to obtain its pivot consistent closure, which is presented in the next section. This filtering 
is achieved w.r.t. a given pivot set P: this is why we now propose the definition of pivot 
consistency with respect to a pivot set P: 

Definition 3.5 Given a CSP V = (X, D,C, R), a root set TZ of (X, C), an TZ-compatible 
assignment ordering x 1, x 2, . . . , x n and a constraint set P C C , we say V is pivot consistent 
with respect to P and this ordering iff P is a pivot set ofV. 

3.2 Pivot Consistent Closure 

Let us assume a CSP V is not pivot consistent, and that we wish to make it so. We then 
obtain a new problem, say V p , but we want this problem to meet some properties. The 
purpose of this section is to specify these properties. We first define the pivot consistent 
closure of a given CSP: 

Definition 3.6 V p = (X,D P ,C P ,R P ) is called the pivot consistent closure of the CSP V = 
(X,D,C,R) iff 

1. V p < V (i.e., V'i E l..n, D p C D i; C C C p , and \/i,j E l..n, ^ C or Cij (£C), and 

2. V p is pivot consistent, and 

3. V p is maximal: there does not exist V' pivot consistent such that V p <V'<V and 

Before presenting two properties of the pivot consistent closure V p of a CSP V, we introduce 
the following lemma, which we will use later in one of the proofs. 

Lemma 3.1 Let V = (X,D,C,R) be a CSP. Let two CSPs V 1 = {X,D l ,C l ,R l ) < V 
and V 2 = (X,D 2 ,C 2 ,R 2 ) < V be pivot consistent w.r.t. a root set TZ, a pivot set P = 
{ x Origin(k) ~^ x k, > \TZ\\ and an TZ-compatible assignment ordering. Let the CSP V s = 
(X, D 3 = D l U D 2 , C 3 = C l n C 2 , R 3 = R 1 U R 2 ), with D 3 = D l U D 2 = {d 3 = d] U d 2 
Vi E l..n} and R 3 = R 1 U R 2 = {r^ = U r\ Vi, j E l..n s.t. c| E C 3 }. Then firstly V 3 
is pivot consistent w.r.t. these same sets TZ and P and this same TZ-compatible assignment 
ordering, and secondly V 1 ,V 2 < V 3 < V . 
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Proof: We first show that V 3 is pivot consistent. Let us denote x 0r igin(k) ~^ x k ( or j 
more briefly, x (k)^>x k ) the pivot whose target is x k in V 1 and V 2 , for all k > \1Z\. Let 
us show that for all k > \1Z\, x Q [k)^x k ls the pivot of Y k -\ m V 3 ■ Let Xj in Yk-\ such 
that c jk G C 1 n C 2 . For all {a j: a 0{k) ) G R 3 j0 (ky either ( a j, a O(k)) £ R )o(k) or ( a ji a O(k)) G 
RjO(k)- ^ ( a ji a o(k)) £ - R )o(A;) : as * s pi v °t consistent, there exists in such that 
{aj,a k ) G (therefore fl 3 ^,) and {a ^,a k ) G (therefore R 3 ^ k)k )- It is the same if 

{a,j,a (k)) G B jO(fc)- c i*: and x O(*0 x ^ are consequently x^-compatible: x ( A .) x fc is 
pivot of Yk-i in V s . 

For all A; > \TZ\, there exists — > Xk pivot of ^ 3 i s pivot consistent. 

Let us now show that V 1 ,V 2 <V 3 < P: 

• V-i G l..n, £)f = U £>■ , so r>- C d 3 and r> 2 C of. (1) 
Furthermore, we have d\ C r>j and £>f C n^; consequently, r> 3 C Dj. (1') 

• C 3 = C 1 n C 2 , so C 3 C C 1 and C 3 C C 2 . (2) 
Also, since C C C 1 and C C C 2 , we have C C C 3 = C 1 n C 2 . (2') 

• V-i, j G l..n such that c| G C 3 , = R-j U so R^ C fj?. and i*?- C fj?.. (3) 

• Vi, j G l..n such that G C, do we have R 3 j C fj^? Let G C. Then cf- G C 3 
(C C C 3 ). By construction of R 3 , we therefore have R 3 j = r\^ U R 2 j. Now fljj C i^- 
and B 2 ^ C which implies fi^- C fl^. (3') 

We deduce from (1), (2) and (3) that V 1 < V 3 and V 2 < V 3 , and V 3 < V from (1 ; ), (2') 
and (3'). □ 

We can now present the properties we announced in the beginning of this section. 

Property 3.2 There only exists one maximal V p < V pivot consistent, for a given root set 
1Z, a given pivot set and a given assignment ordering. 

Proof: Suppose that V p is not unique: thus, there exists V = (X, D' , C" , R') < V such 
that V is maximal and pivot consistent, and V' / V p . Let us build the CSP V" = (X, D p U 
D',C P n C\R p U R'). This CSP is pivot consistent and V p < V" < V (Lemma 3.1): V p is 
not maximal, which contradicts the previous assumption. □ 

The second property of V p is perhaps the most important: 

Property 3.3 Both CSPs V and V p are equivalent, in that they have the same set of solu- 
tions. 

Proof: Let us denote S(V) and S(V P ) their respective sets of solutions. Obviously, S(V P ) C 
S{V). Let us now show that S{V) C S{V P ): suppose there exists I n = (si, . . . , s n ) solution 
of V which is not a solution of V p . Then G l..n such that (si,Sj) R?j. Now 

(si,Sj) G Rij (I n is a solution of V): as "P^ is maximal, that means that (si,Sj) has no 
support for at least one constraint: 3k s.t. Va^ G (si,a k ) R%k or (sj,a,k) flj/i- We 
therefore have (sj, s k ) fiifc or (sj, s k ) iij^: 7 n is not a solution of V, that is S(P) C S(V P ). 
□ 
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We now know that there exists a unique CSP V' , called the pivot consistent closure of 
V, such that V < T 7 , T 7 ' is maximal and T 7 ' is pivot consistent. Moreover, this CSP has the 
same set of solutions as V, the problem it comes from. But, even though it is interesting to 
know that such a CSP exists, it seems reasonable to wish to obtain it: that is the purpose 
of the next section, which proposes an algorithm achieving a pivot consistent filtering. 

3.3 A Filtering Algorithm 

First, a notation: the set of the functional constraints x^ — > x& chosen to be the pivots of 
each Yk is called a set of pivot candidates, and denoted PC. After the filtering, PC = P, the 
pivot set. We will suppose in the course of this section that both the assignment ordering 
and PC are known; we explain how they are obtained in appendix A. 2. 

The algorithm is composed of several procedures. Their different levels represent the 
three steps followed to define pivot consistency. We present these procedures in an ascending 
way: 

1. Make the constraints Cjk and c^k = X/j — > x& x^-compatible. 

2. Compute a pivot — > x^ for all x^ in X — 1Z 

3. Achieve pivot consistency for the CSP. 

Procedure Compatible (x^, x x j) makes constraints x^ — > x^ and Cjk x^-compatible: it 
removes from R^j those tuples which do not have a common support in for these two 
constraints. If necessary, it creates the constraint C'hj between the variables x^ and Xj. 

Procedure Compatible (x/j, x k , Xj) 
begin 

for all (ah,aj) € Rhj do 
at <~ fhk(ah) Q 

if (a h ,a k ) & Rhk or (aj,a k ) g R jk then Suppress (a h ,aj) from R hj 
end for 
end 

Procedure Pivot (x/^x^) makes all constraints Cjk containing x& and such that Xj 6 Yjfc-i 
Xfc-compatible with X/j — > x^ by successive calls to the subroutine Compatible (x/ t , X&, Xj). 
After its computation, x^ — > x^ is therefore a pivot of i^-i- 

Procedure Pivot (x/,, x k ) 
begin 

for all Xj (E Yk-\ s.t. Cjk G C do Compatible (x/j, x k , Xj) 
end 

The CSP is pivot consistent iff any variable x^ of X — 1Z is the target of a pivot of l^-i- 



9. By convention, if ah has no support in Du for Chk, we denote fhk(dh) = e, and, of course (a,h,s) $ Rhk 
and (a,-, e) £ Rjk- 
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Pivot consistent filtering 
begin 

for k <— n to r + 1 do 

Xh i— Origin of the pivot candidate whose target is Xj. 
Pivot (x h ,x k ) 
end for 
end 

Proposition This algorithm achieves pivot consistency (the resulting problem is the pivot 
consistent closure of V w.r.t. the chosen data). 

Proof: We clearly see that after computing Pivot (x/^x^), — > x& is a pivot of Y^-i- 
Moreover, only the necessary suppressions are performed: according to the definition of 
Compatible, a pair is suppressed from a relation only if it has no support in the target of 
the current pivot candidate. 

We consequently simply have to show that any suppression due to Pivot (x/^x^) has 
no influence on the pivots x g — > X\ for I > k. More precisely, if x g — > X\ is a pivot of Y\-\ 
before computing Pivot (x/^x^), then it will remain so. 

For the sake of brevity, we will denote in the course of this proof Pv(k) the computation 
of the procedure Pivot (x/^x^). Let I be such that x g — > X[ is a pivot of Y;_i before 
Pv(k). Suppose x g — > X{ is no longer a pivot of Y\-\ after Pv(k). Consequently, there 
exists j < I such that Cji and c g i are not x^-compatible: B(aj,a g ) G Rj g s.t. (aj,ai) Rji 
or (a g ,ai) R g i. Pv(k) has thus suppressed from Rji or (a g ,ai) from R g i. This is 

impossible, since by definition procedure Pivot (x^x^) only modifies constraints C^, with 
h,i < k, so a fortiori h,i < I (k < I). So, Pv(k) does not influence pivots x g — > X[ for 
/ > k. This algorithm therefore computes the pivot consistent closure of a CSP. □ 

Notice that unlike arc or path consistencies which only need knowledge about the CSP, 
the pivot consistent closure is in addition defined w.r.t. a root set, a pivot set and an 
assignment ordering: Figure 7 shows how the choice of the pivot (s) may influence the 
filtered problem. 

3.4 Application to the Example 

In this section, we illustrate the filtering algorithm using the travel agency example. We 
first recall the initial problem and expose the data necessary to the algorithm (root set 1Z, 
set of pivot candidates PC and 7^-compatible ordering). 

Tl = {GUIDES, CITIES} 

PC = {CITIES COUNTRIES, 

COUNTRIES CURRENCIES, 
COUNTRIES LANGUAGES} 

An 7^-compatible ordering is: 

GUIDES, CITIES, COUNTRIES, CURRENCIES, LANGUAGES 

In order to achieve pivot consistency, we need to compute Pivot (COUNTRIES, LANGUAGES), 
Pivot (COUNTRIES, CURRENCIES) and Pivot (CITIES, COUNTRIES). 
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K = {x 1 ,x 2 ,x 3 } 

The three possible pivots are X\ — > X4, x 2 — > X4 and X3 — > X4. 
The assignment ordering we have chosen is X\, x 2 , x$, x 4 . 

If we choose the pivot x 2 — > X4 the filter- 




Figure 7: Pivot consistency: influence of the pivot 



CURRENCIES LANGUAGES 




Figure 8: Constraint and consistency graphs of the initial V, 
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Figure 9: After Pivot (COUNTRIES, LANGUAGES) 



We now detail these computations: 

• Computation of Pivot (COUNTRIES, LANGUAGES) (Figure 9): 

We perform Compatible ( COUNTRIES, LANGUAGES, COUNTRIES) which does not 
modify anything, and Compatible (COUNTRIES,LANGUAGES, GUIDES), which cre- 
ates the constraint {GUIDES, COUNTRIES}. 

• Computation of Pivot (COUNTRIES, CURRENCIES): 

We perform Compatible (COUNTRIES, CURRENCIES, COUNTRIES) which does not 
modify anything. 

• Computation of Pivot (CITIES, COUNTRIES) (Figure 10): 

We perform Compatible (CITIES, COUNTRIES, CITIES) which does not change any- 
thing, and Compatible (CITIES, COUNTRIES, GUIDES), which modifies the constraint 
{GUIDES, CITIES}. 

3.5 Pivot Consistency vs Path Consistency 

We previously said that arc-and-path consistency was not necessary to give some properties 
to a functional CSP, and that pivot consistency was sufficient. We first show in this section 
that pivot consistency is a weakened form of path consistency, and in section 4 we will 
present some properties and a method for solving functional CSPs based upon them. 

Note that, since pivot consistency is directed, one could wish to compare it not only with 
"full" path consistency, but also with directional path consistency. We therefore propose to 
present relationships between pivot consistency and both path consistency and directional 
path consistency. We will thus see that pivot consistency is a restricted version of directional 
path consistency as well. 
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Figure 10: After Pivot (CITIES, COUNTRIES): end of the filtering 



First, let us recall a remark we made when defining x^-compatibility: 

A CSP is path consistent iff for all x& in X and for all X/j,Xj in X, {x/j,x&} 
and {Xj,Xfc} are x&- compatible. 

A slight difference appears with directional path consistency: 

A CSP is directional path consistent iff for all x^ in X and for all X/j,Xj in X 
such that h,i < k, {x/j,x&} and {Xj,x&} are x^- compatible. 

Pivot consistency does not need so many requirements. First, ^-compatibility is only 
required for variables in X — 1Z. Furthermore, for each of these x&, not all of the possible 
pairs of constraints are filtered: as this consistency is directed, x^-compatibility is only 
achieved for constraints whose first variable precedes x^ (the second one) in the assignment 
ordering. Finally, the filtering is only achieved w.r.t. the pivot candidate whose target is 
Xk (say Xh k — > Xk). In other words, 

A CSP is pivot consistent iff for all x^ in X — 1Z and for all Xj in X such that 
i < k, {Xh k , x&} and {Xj, x&} are x&- compatible, provided that Xh k — > x& 6 PC. 

To sum up, one can roughly say that for each x&, path consistency needs to consider 0(n 2 ) 
triangles, directional path consistency 0(k 2 ) and pivot consistency O(k). The difference 
between path and pivot consistencies clearly appears on the travel agency example: only 
two constraints are altered. The first one ({GUIDES, COUNTRIES}) is created, and the 
second one ({GUIDES, CITIES}) is modified, whereas path consistency modified the five 
existing constraints, and created the five other possible constraints, making the constraint 
graph complete (see Figure 11). However, this example does not highlight any difference 
between pivot consistency and directional path consistency. Figure 12 presents how pivot 
consistency and directional path consistency may differently modify the constraint graph 
of a CSP when processing a node x&. 
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Figure 1 1 : Path consistent closure of V e 
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x g Directional path consistency 





Figure 12: Processing pivot consistency and directional path consistency on node x^, (with 
g, hk,i, j < k): the bold edges represent the constraints possibly created by pivot 
consistency (w.r.t. x^ k — > x^), the dotted edges the extra ones possibly created 
by directional path consistency 
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Time complexity is lower for pivot consistency than for path consistency or directional 
path consistency, as will now be shown. 

3.6 Complexity 

Procedure Pivot is computed exactly once for each variable x k in X — 1Z; it then calls 
at most k — 1 times procedure Compatible (once for each of the k — 1 former variables). 
Achieving x^-compatibility between two constraints C'h k k and C\ k corresponds to achieving 
path consistency for relation w.r.t. variable x k . This operation needs 0(d 3 ) in general. 
But since the constraint Ch kk is functional, this complexity is now 0(d 2 ). The overall 
complexity of the filtering is therefore 

n 

0{ ( k - !M 2 ) = 0((n 2 - r 2 )d 2 ) 

(instead of 0(n 3 d 3 ) for path consistency or directional path consistency 10 ; n is the number 
of variables, d is the size of the domains and r is the size of the root set). 

4. From a Local Instantiation to a Solution 

The intrinsic characteristic of a local consistency is to ensure that a partial instantiation 
can be extended to a new variable. We first present the properties of pivot consistency, in 
particular the conditions under which a consistent partial instantiation may be extended to 
a solution. We then explain how to compute the data required by the filtering algorithm, 
and finally present a method for solving functional CSPs. 

4.1 Properties 

We proceed in two stages: first, the addition of a new variable to the current instantiation, 
then the extension from the root set to a solution. 

Property 4.1 Let V = (X, D,C, R) be a CSP, I k -\ be a consistent instantiation of Y k _i, 
and x k be a new variable to instantiate. If there exists x^ 6 Y k _\ such that x^ — >■ x k is a 
pivot ofYk-\, then I k -\ may be extended into I k , a consistent instantiation ofY k . 

Proof: Let us show that any constraint included in Y k is satisfied. Let x^ 6 Y^-i such 
that Xh — > Xk is a pivot of and its value in Ik-i- Consequently, there exists 

a k = f h ^ k (a h ). Let us denote I k = (a u ... ,a h ,.. .,a k -i,a k ). 

1. Any constraint satisfied by I k -\ and which does not contain x k is therefore obviously 
also satisfied by I k . 

2. For all i < k s.t. c\ k 6 C. First, Xj 6 Y k _\: let ai be its value in I k -\. Second, 
Xh — > x k is a pivot of Yjt-i- c ik an d x^ — > x k are consequently v^-compatible. We 
thus have (ai,a k ) 6 Ri k ; so, c\ k is satisfied by I k . 

Every constraint included in Y k is satisfied by I k . Hence I k is a consistent instantiation. □ 

10. Note that achieving directional path consistency on functional CSPs is still 0(n 3 d 3 ) and not 0(n 3 d 2 ), 
since first functional constraints are not necessarily directed the right way (where a functional constraint 
is "properly" directed if its origin is before its target in the instantiation ordering), and second non 
functional constraints may be created by the processing. 



465 



David 



We now present the theorem at the center of the method we introduce further. 

Theorem 1 Let V = (X, D,C, R), TZ C X a root set, and an 1Z- compatible assignment 
ordering such that V is pivot consistent w.r.t. them. If a consistent instantiation of TZ 
exists, then it can be extended to a backtrack-free solution. 

Proof: The variables of X — 1Z only need to be instantiated along with the 7£-compatible 
ordering: as each of them is the target of a pivot (since V is pivot consistent), Property 4.1 
is verified step by step from the instantiation of TZ to a solution. □ 

4.2 A Solving-by-Decomposing Method 

In section 4.1 (Theorem 1), we saw that, given a consistent instantiation of TZ, and provided 
that the problem is pivot consistent w.r.t. an 7^-compatible assignment ordering and a pivot 
set, we can not only guarantee there is a solution to the whole problem, but also find it 
without any backtracking. We can therefore deduce a method from this property. We now 
introduce this method, decomposed into four phases, and for each of them, we present the 
time complexity and the result of its computation on the example. 

Phase 1 Computation of a root set TZ, a pivot candidate set PC and an 7£-compatible 
assignment ordering. 

This phase can roughly be outlined as finding the sources of a graph where each 
strongly connected component is reduced to one node (computation of TZ), finding a 
set of directed arcs covering all nodes in X — TZ (PC) and then performing topological 
sort (assignment ordering). For more details, see appendices A.l and A. 2. 

Complexity: 0(e + n) to compute TZ (appendix A.l), and also 0(e + n) for PC and 
the assignment ordering (appendix A. 2). So, for the whole phase: 0(e + n). 

The root set is TZ = {GUIDES, CITIES}. 

The pivot candidate set is PC = {CITIES COUNTRIES, COUNTRIES CURREN- 
CIES, COUNTRIES LANGUAGES} 

The 7£-compatible assignment ordering we chose is 

GUIDES, CITIES, COUNTRIES, CURRENCIES, LANGUAGES. 

Phase 2 Pivot consistent filtering 

Complexity: 0((n 2 — r 2 )d 2 )) (see 3.6) 

The pivot consistent problem is computed in section 3.4. 

Phase 3 Instantiation of TZ 

Complexity: 0(en,d r ), where en is the number of constraints included in TZ. 

There are five consistent instantiations of TZ = {GUIDES, CITIES}: 

{(Alice, Paris), (Bob, Paris), (Bob, London), (Bob, Washington), (Bob, New- York)} 
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Phase 4 Instantiation of X — TZ 

Complexity: 0(n — r): extending a consistent instantiation of TZ to a solution only 
needs linear time. 

The solutions of the problem are: 

(Alice, Paris, France, FrF, French) 
(Bob, Paris, France, FrF, French ) 
(Bob, London, GBob,£, English) 
(Bob, Washington, USA, $, English) 
(Bob, New-York, US A, $, English) 

The overall time complexity of this method is therefore 0((n 2 — r 2 )d 2 + e-]id r ). This method 
thus reduces the search for one solution (or all) of a functional CSP to the one induced by 
its root set TZ. As the time required obviously depends on the size of TZ, we immediately 
see how interesting it is to compute a minimum sized root set. A direct consequence is that 
the smaller 7Z is, the more efficient this method will be. It consequently seems worthwhile 
to know its size as soon as possible: this has be done in phase 1. 

A last remark before concluding: the reader probably noticed that the number of con- 
sistent instantiations of TZ is the same as the number of solutions of the whole problem. 
This is not a coincidence, as we see below: 

Property 4.2 Let V be pivot consistent, and let TZ be its root set. The number of solutions 
of V equals the number of consistent instantiations of TZ: each solution is the extension of 
exactly one of those instantiations. 

Proof: Let I r = (ai, . . . , a T ) be a consistent instantiation of the root set TZ = {x\, . . . ,x r }. 
According to Theorem 1, can be extended to at least one solution to V. TZ is a root 
set of Q. Any in X — TZ is therefore a descendant from a variable x Tk of TZ. There 
consequently exists a unique value 6 such that {a rk ,a,k) 6 R rk k-> and, so, an only 
consistent instantiation of TZ U x& including I r , therefore only one solution including I r . 
Moreover, as two different instantiations of TZ obviously cannot be extended to the same 
solution, there consequently exists a one to one mapping between the set of the consistent 
instantiations of TZ and the set S of the solutions to V. □ 

5. Conclusion 

Taking into account semantic properties of the constraints is a still recent approach to 
increase the efficiency of finding solutions to CSPs. This paper belongs to this frame. 

We first introduced a new local consistency, pivot consistency, to deal with functional 
constraints for a lower cost than path consistency (we furthermore showed that pivot con- 
sistency is a weak form of path consistency). We then proposed an algorithm to achieve 
pivot consistent filtering. Later, we studied some properties that pivot consistency provides 
a functional CSP with, in particular conditions under which a consistent instantiation can 
be extended to a solution. 

We were then led to present a decomposition method based on those properties, which 
decreases the complexity of solving functional CSPs. An interesting point is that the search 
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for solutions to the former problem is reduced to the search for solutions to the subproblem 
induced by a particular subset, the root set; this new problem can be solved by any method, 
including heuristics. 

Furthermore, we must add that this method deals with problems in which not all the 
constraints have to be functional (as we saw in the example). This is in our opinion an 
interesting improvement: most previous work dealing with properties of specific classes of 
constraints assumes that all the constraint of the network possess a given property, and 
only applies in this case. 

However, some problems remain. First, and this is a classical problem as far as pre- 
processings are concerned: in what extent is the application of the method we present here 
useful? A partial answer can be given to this question: we indeed saw that one of the 
advantages of this method is that the cost of the search is known early in the process; one 
can thus estimate if the continuation of the process is worthwhile or not. 

On the other hand, pivot consistency and the associated consistency properties can of 
course be generalized to non binary CSPs; yet, some problems appear: first of all, finding a 
minimum root set becomes exponential (it is the same problem as finding a key of minimum 
size in the field of relational databases, Lucchesi Sz Osborn, 1978); moreover, as we have 
seen, some constraints may be induced by the filtering processing; if this is not a problem 
for binary CSPs (the new constraints are still binary), in n-ary CSPs, some constraints may 
induce new ones of greater arity, which may from step to step lead to an explosion of the 
CSP arity. 

Finding and characterizing classes of problems having root sets of small size as well as 
problems whose arity does not grow when achieving pivot consistency are therefore in the 
continuity of the work we presented in this paper. 

To conclude, taking into account the semantics of constraints seems to be an interesting 
research area, and, as such, we can think of extending this kind of study (characterization of 
properties and of processings specific to a class of constraints) to other classes of constraints. 

Appendix A. Computation of Structural Data 
A.l Computation of a Minimum Root Set 1Z 

The example used in this paper is not suited to the illustration of the computation we 
present here; we therefore use another one (Figures 13 to 15). 

Computation of Q q , reduced graph of Qf. 

1: Compute T = {t\, . . . , T p }, the set of the strongly connected components of Qf. 

2: Compute graph Q q = (T,C q ), where each vertex ti of T is the reduction of a 
strongly connected component Tj, and there exists an arc from ti to tj if and 
only if an arc existed from a vertex of Tj to a vertex of Tj . 

Computation of the sources of Q q , and the root set of Qf. 

3: Compute the set S q = {t Sl , . . . , t Sr } of the sources of Q q . 

4: Choose for each source t Sk of S q any vertex of T Sk , which will represent it. 

5: The set 1Z processed in that way is a root set of Qf. 
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Tl = {X L } 

T 2 = {X 2 ,X 3 ,X 4 } 

T 3 = {X 5 } 

Ti = {X 6 } 

T 5 = {X 7 ,X 8 ,X 9 } 



Figure 13: Step 1: the strongly connected components 




Sq — {^l,^} 



Figure 14: Steps 2 and 3: the reduced graph and its sources 




The only possible vertex to repre- 
sent t\ is X\. 

A vertex is chosen amongst 
T2 = {^2,^3,^4} to represent ti, 
for instance Xi- 

We thus obtain 1Z = {x\, x 2 } 



Figure 15: Steps 4 and 5: a minimum root set 
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A question arises: 7s the set TZ really a minimum root set? 

There are actually two questions: first, is TZ a root set? and if so, is this set minimum? 

We first show that the set TZ we computed is a root set: 
Let us show that any vertex is a descendant from an element of TZ. Let Xj be a vertex of 
X — TZ. Two possibilities: 

1. An element of 1Z belongs to the same strongly connected component as Xj. By defi- 
nition, Xi is then a descendant from this element. 

2. No element in 1Z belongs to the same strongly connected component as Xj. Let us 
denote Tj the strongly connected component which contains Xi, and ti its reduction in 
the graph Q q . So, ti is not a source of this graph; consequently, there exists a source 
t Si in Q q such that ti is a descendant from t Si . Let x Si be the vertex of X we chose to 
represent t Si in TZ. According to the definition of the reduced graph, knowing that ti 
is a descendant from t Si , any element of Tj is a descendant from any element of T Si , 
and therefore from x Si . Xj is consequently a descendant from an element of TZ. 

TZ is a root set. □ 

We now show that TZ is minimum (i.e., there does not exist any root set of smaller size): 
Let us denote r = \TZ\. Assume there exists a root set TZ 1 whose size is r' < r. Let SCC(TZ') 
be the set of the Strongly Connected Components that contain all the elements oiTZ', and 
let T(TZ') be the set of their reductions in Q q . By definition of TZ, the reduced graph Q q has 
r sources. There exists at least one source t Sk in Q q which does not belong to T(TZ') (since 
|7~(7£')| <r'<r). Consequently, there exists no path from an element of T(TZ') to t Sk . 
Let us denote SCC(k) the strongly connected component reduced to t Sk . Prom the definition 
of the reduced graph, there exists no path from an element of TZ 1 to an element of SCC(k). 
Consequently, TZ' is not a root set. 

TZ is thus a minimum root set. □ 

Complexity: It is the same as the complexity required for computing the strongly connected 
components, that is 0(e + n) (Tarjan, 1972), if e is the number of edges and n the number 
of vertices. 

A. 2 Choice of the Pivots and Computation of an 7^-compatible Order 

We now describe how to choose the pivot candidates and to compute an 7£-compatible 
ordering. We first present the conditions the pivot candidates must satisfy; we then present 
an algorithm that computes both PC and an 7^-compatible ordering. The definition of 
pivot consistency implies two conditions on the pivots: 

1. Any Xfc in X — TZ must be the target of one and only one pivot, and no variable of TZ 
is the target of a pivot. 

2. If X/j — > Xfc is a pivot, then X/j is before x& in the ordering. 

For each X& in X — TZ, we have to choose a variable X/j in X so that X/j — > x& (x/j necessarily 
exists, since x^ does not belong to TZ). Moreover, PC must not contain any circuit, which 
would be in contradiction with condition 2. A way to prevent circuits is to mark every 



470 



Pivot Consistency for Functional CSPs 



selected variable, and to select a new variable among the unmarked ones, which is the 
target of a functional constraint whose origin is already marked. This functional constraint 
is then included into PC, and its target (the new variable) is then marked. The set PC we 
obtain this way is the set of the pivot candidates used by the pivot consistency algorithm 
(section 3.3). This set induces a partial order -<n on X (we only have to add the transitivity 
constraints). In order to satisfy the conditions of Definition 3.3, the assignment ordering 
consequently has to be a linear extension of the partial order -<n prefixed by TZ. Actually, 
computing -<n is not necessary: PC is sufficient to compute the linear extension. There are 
no conditions on the variables of 1Z; the algorithm can thus be decomposed into two steps: 

1 . Ordering of TZ 

Number from 1 to r = \TZ\ the variables of TZ and mark them 

2. Ordering of X — 1Z 
Repeat 

Choose a new unmarked variable in X — 1Z s.t. there exists a marked x^ s.t. 

Number and mark x & 
Until all variables are marked 

Figure 16 presents an algorithm computing both the set of pivot candidates PC and an 
"/^.-compatible assignment ordering. Let us repeat that any linear extension issued from PC 
is an 7£-compatible ordering. 

• The set Marked contains the variables that have already been treated 

• The set NextPossible contains the next variables that can be chosen, that is the 
unmarked ones which are the target of a functional constraint whose origin is marked 

• The sets Origin [j] contain the origins of functional constraints whose targets are 
described above (NextPossible) 

• Num represents the number of the current variable in the assignment ordering 

Correctness: 

• Set PC: 

Condition 1: a functional constraint — > x^ is added to PC each time a new variable 
Xk is chosen: we therefore need to prove that any variable of X — TZ is selected exactly 
once, and that none is selected from 1Z. The variables are selected from NextPossible, 
which only contains unmarked variables; since no marked variable may be unmarked 
back, and every selected variable is marked, every variable can be selected at most 
once. 

Let us now prove that every variable of X — TZ is actually (at least once) selected: at 
any time, NextPossible is the set of all direct descendants of the marked variables 
that do not belong to TZ; at each outer loop of step 2, a new variable is extracted 
from NextPossible and marked; from the definition of TZ, all variables of X — TZ will 
therefore be reached and inserted in NextPossible. 
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begin 

Initializations 

Marked <- 
NextPossible -5— 

for each i in l..n do Origin[i] •<— 
Num <- 1 

Step 1: Numbering of TZ 

while Marked ^ TZ do 

Choose Xj unMarked in 7?. 
Mark Xj (* Add Xj to Marked *) 
Number Xj with Num 
Add 1 to Num 

for each Xj — > Xj s.t. Xj is unMarked and Xj ^TZ do 

Add Xj to NextPossible 

Add Xj to Origin [j] 
end for 
end while 

Step 2: Numbering of X — TZ — Construction of PC 
PC ^ 

while Marked / I do 

Choose Xfc in (and remove it from) NextPossible 
Mark x& (* Add x& to Marked *) 
Number x& with Num 
Add 1 to Num 

Choose X/, in (and remove it from) Origin [A,] 
Add X/j — > Xk to PC 

for each x^— > x p s.t. x p is unMarked do 

Add x p to NextPossible 

Add Xfc to Origin [p] 
end for 
end while 

end 

Figure 16: Computation of a set of pivot candidates PC and an 7^-compatible ordering 

Condition 2: by construction of the sets Origin [j] , a variable Xj is inserted in 
Origin [j] immediately after Xj has been numbered; this only applies for the un- 
marked Xj. So, Xj is necessarily numbered after Xj. 

• Is this assignment ordering 7^-compatible? 

The first variables are obviously those taken from TZ (step 1). The other condition is 
a direct consequence of condition 2 above. 
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Complexity: 

Both steps have the same structure: choose a new variable Xi, mark and number it, and 
for each of its (unmarked) direct descendants, do some 0(1) operations (during step 2, 
a constraint is added to PC which also requires 0(1)). The overall time complexity is 
therefore 

n 

0(£(l + |r + (t)|)) = 0(n + e) 

2 = 1 

where T + (i) is the set of the direct descendants from Xi. 
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